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The physics involved in the fundamental conservation equations of the spin and orbital angular 
momenta leads to new laws and phenomena that are disclosed here. To this end, we analyse the 
scattering of an electromagnetic wavefield by the canonical system constituted by a small particle, 
which is assumed dipolar in the wide sense. Specifically, under quite general conditions these 
laws lead to understanding the contribution and weight of each of those angular momenta to the 
electromagnetic torque exerted by the field on the object, which is shown to consist of an extinction 
and a scattering, or recoil, part. This leads to an interpretation of its effect different to that taken 
up till now by many theoretical and experimental works, and implies that a part of the recoil 
torque cancels the usually called intrinsic torque which was often considered responsible of the 
particle spinning. In addition, we obtain the contribution of the spatial structure of the wave to this 
torque, unknown to this date, showing its effect in the orbiting of the object, and demonstrating 
that it often leads to a negative torque on a single particle, i.e. opposite to the incident helicity, 
producing an orbital motion contrary to its spinning. Furthermore, we establish a decomposition of 
the electromagnetic torque into conservative and non-conservative components in which the helicity 
and the spin angular momentum play a role analogous to the energy and its flux for electromagnetic 
forces. These phenomena are illustrated with examples of beams, also showing the difficulties of 
some paraxial formulations whose fields do not hold the transversality condition. 

PACS numbers: 33.57.-l-c, 42.25.Ja,75.85.+t 


INTRODUCTION 

There is a growing interest in phenomena related to 
the angular momenta of light and other electromagnetic 
fields. Advances in particle manipulation are making ex¬ 
periments, their assessment, interpretation with theories, 
and applications of increasing study and capability. This 
has opened a new area of research related to the twisting 
of both the polarization and the wavefronts of electomag- 
netic fields |ll42Cll| , and their effects on matter 2l|,[2^ [24- 

lillilil 

On the one hand, conservation laws for the helicity 
and spin of wavefields were recently established 

based on initial extensions [llj| of known conservation 
equations of electromagnetism. Such laws appear as fun¬ 
damental as those for the energy, and linear and angular 
momenta. In particular Q put forward the conservation 
of spin and orbital angular momenta, separately. On the 
other hand, studies on the momenta transfer to matter, 
and their mechanical action on objects, are enlarging the 
area of optical manipulation of objects. Among these, 
new bizarre effects like negative optical torques (i.e. op¬ 
posite to the helicity of the illumination) on sets of parti¬ 
cles have been predicted [H, [s^l with circularly polarized 
plane waves, as well experimentally observed with Gaus¬ 
sian beams on extended objects 3^. This phenomenon 
keeps an analogy with pulling forces 3^ U | that so much 
attention have recently attracted. 

However, theory and experiments on electromagnetic 
torques are less developed and understood than their 


counterparts of optical and electromagnetic forces [3S 
54| . Observations are more difficult to control and to 


quantitatively interpret with existing models. With few 


exceptions [34 [^ , most experimental [27|, [2g and theo¬ 
retical [ 29 I - I 34 1531 studies employ a static formulation, 
(perhaps following the path of pioneering work 2l|), 
which was shown 37] to be incomplete and not com¬ 
patible with energy and angular momentum conserva¬ 
tion. Only for extremely small (i.e. Rayleigh) particles 
in terms of the wavelength the static approximation is 
valid; this is further discussed in Section 10, where the 
total electromagnetic torque is addressed. 

Although several reports 


44 ] realize 


the need of energy absorption by the object in order that 
it experiences a torque, no explicit demonstration exists 
of the role played in this effect by the variation of incident 
spin and orbital angular momenta, even thou^ they are 
calculated in some cases after scattering [1,H, 0) 0 - 
Moreover, studies based on the static approximation |27l - 
0) 32, only deal with the so-called intrinsic torque 


which, as as shown below cannot account for the angular 
momentum transfer nor can describe the resulting torque 
experienced by the object through energy absorption. 


In a recent work we established the significance 
of the conservation of electromagnetic helicity for scat¬ 
tering objects and its relevance for energy transfer be¬ 
tween small particles and molecules, as well as for circular 
dichroism. Now the aim of this paper is two-fold: On the 
one hand, we integrate the conservation laws for the spin 
and orbital angular momenta, extracting their physical 
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meaning for the interaction of fields and objects, specifi¬ 
cally in the scattering of a wavefield by a small particle. 
On the other hand, we predict the contributions of each 
of these angular momenta conservation laws to the opti¬ 
cal torque exerted by the field on the object; also taking 
into account that part due to the spatial structure of the 
field. This is illustrated with a generally magnetodielec¬ 
tric bi-isotropic dipolar particle in the wide sense, 

i.e. whose scattering is fully described by the first order 
partial waves, (namely, Mie coefficients if it were a sphere 
[54ISTI . l58|l. This analysis opens a collection of new phe¬ 
nomena related to the electromagnetic torque, suceptible 
of further investigation from the theoretical and experi¬ 
mental points of view. Our discussion is classical, how¬ 
ever if one considers the quantum nature of light, these 
results may be obtained in terms of the number operators 
associated to left and right circular polarized modes, ex¬ 
pressed by the product of the corresponding annihilation 
and creation photon operators, (see |59|-l6lj), specially 
since methods to identify twisted beams with different 
topological charge have been devised [Hll - I^ . 

A Coulomb gauge is chosen dealing with fields that sat¬ 
isfy the transversality condition. This allows a separation 
of the spin and the orbital angular momenta, thus confer- 
ing sense to the existence of their respective conservation 
laws [3,[i3, separately. It is known, however, that cer¬ 
tain formulations of optical beams have a problem with 
the transversality condition, and this is also noticed in 
this work. Hence, within this context, and within the va¬ 
lidity of these conservation laws, we reach the following 
results: 

1. The conservation laws for the spin and orbital an¬ 
gular momenta lead to their respective contributions to 
the time-averaged electromagnetic torque on the object, 
which from each of them is composed of two parts: (a) 
A extinction torque due to the extraction of the corre¬ 
sponding (spin or orbital) angular momentum from the 
incident wave by scattering from the body, described by 
the interference between the incident and the scattered 
fields, (b) A scattering or recoil torque due to the scat¬ 
tered field interfering with itself. 

2. The recoil torque due to the orbital angular mo¬ 
mentum conservation is equal to that yielded by the spin 
angular momentum conservation; both of them add to 
produce the total recoil torque, which is the same as 
that derived from the conservation of total angular mo¬ 
mentum , and coincides with its emission rate by the 
electric and/or magnetic dipole induced on the particle, 
as should be. 

3. The extinction torque obtained from the spin an¬ 
gular momentum conservation contains two terms: one 
is the so-called intrinsic torque by several authors so 
far, usually employed in most studies on dipolar parti¬ 
cles based on the static theory, and which however we 
show to be cancelled by a part of the total recoil torque, 
which is also intrinsic as it does not depend on the ori¬ 


gin of coordinates; the other term comes from the spatial 
structure of the incident field which, in turn, is equal 
to twice the corresponding extinction torque yielded by 
the conservation of orbital angular momentum. If the 
incident wave is plane, obviously there is no term due 
to the field spatial structure; hence the orbital angular 
momentum conservation yields no extinction part, which 
is consistent with the absence of incident orbital angular 
momentum in such a wave. Also, in this latter case, the 
above mentioned cancellation conveys a transfer of spin 
angular momentum to the particle accounted for by the 
remaining part of the recoil torque through absorption of 
the incident energy. 

4. Therefore, if the incident field is a plane wave, ac¬ 
cording to Point 3 above the existence of a recoil torque 
stemming from the conservation of orbital angular mo¬ 
mentum has to come from the incident spin angular mo¬ 
mentum, this being a manifestation of the spin-orbit in¬ 
teraction which is included in the conservation laws. 

5. For a circularly polarized incident plane wave, if the 
particle produces a circularly polarized scattered field, 
(for which it is sufficient that its electric and magnetic po¬ 
larizabilities coincide), then the electromagnetic torque 
and force on that particle become proportional to each 
other. 

6 . Taking into account the spatial structure of the 
field, the torques derived from the spin and orbital con¬ 
servation laws admit a decomposition into conservative 
and non-conservative components, analogous to that of 
electromagnetic forces. There is a gradient component of 
the torque, where now the wavefield helicity plays a role 
analogous to that of the energy for the optical force; and 
hence may be employed as a rotational optical tweezer 
that aligns the torque acting on the particle along an 
equilibrium direction. Likewise, there is a component 
given by the spin angular momentum, (which coincides 
with the helicity flow in the dual-symmetric formulation 
for quasi-monochromatic fields Q here employed), analo¬ 
gous to the radiation pressure, or Poynting vector part, of 
the force. In addition, there are other non-conservative 
components that describe torques due to circulation of 
field vortices around the dipolar object. 

7. In contrast with plane waves, incident fields with 
certain spatial structures, and some beams in particular, 
may easily produce negative electromagnetic torques on 
one single small particle. 


THE FLOW OF HELICITY 


We shall address time-harmonic fields with electric 
and magnetic vectors: £{r,t) = Ji[E(r) exp(—fwt)] and 
B{r,t) = 3?[B(r) exp(—fwt)], a functional form that also 
applies to all potentials and currents. 5R denotes real 
part. Although such fields are not source-free 65, 
they will be considered in regions with no sources so that 
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they admit an angular spectrum representation of plane 
wave components (^. 1^. either convergent, divergent, 
or both, satisfying at /cr —>■ oo Sommerfeld’s radiation 
condition. Such fields are of course transversal, namely, 
they hold Maxwell’s equations: V ■ £ — ■ B — Q. 

In a non-absorbing dielectric medium of refractive in¬ 
dex n = y/ejl, (e and /i representing the dielectric per¬ 
mittivity and the magnetic permeability), the helicity 
density Jff, the flow of helicity density and the ten¬ 
sor Mij of density of flow of helicity flux, of this light 
field are defined in the dual-symmetric formulation as 
(cf. 


the transversality condition is not taken into account, 
dealing with some representations of optical beams which 
do not fulfill it, thus preventing their inclusion in such 
decomposition of the angular momentum and, in general, 
in any theory based on such transversality equations. 

From the above conditions, the definitions, ^ and 
jVij hold the continuity equations 0, [I^ 

jT-t V-^ =(5) 

# -f V • = -Si, (6) 


J^ = ^{^A-B-eC-£), ( 1 ) 

^=^{£xA + BxC), ( 2 ) 

1 

Mij = £jCi ~ S\.jBi) + 

s (J-A-B-S-C)], (^,J = 1,2,3). (3) 

en 

In ([1]) and ([2|) A and C are vector potentials such that 
B = V X A and £ = —V x C and are transversal in a 
Coulomb’s gauge [l^: V-M = V-C = 0. From Maxwell’s 
equations: 

• • r Att 

A = -c£, C = - \/xA+—JC. (4) 

e/i e 

The upper dot meaning dt, c being the light speed in 
vacuum. In writting C in (j?]) we have taken into accent 
that the electric current density, which we denote as Sf, is 
transversal since the existence of A and the law V ■ e£ = 
Airp convey that the electric char ge density p and any 
(static) scalar potential are zero [68|, and thus J has 
been expressed as J = V x 1C. 

The transversality of A and C also amounts through 
the two first Eqs. 0 to that of £ and B, and hence is 
compatible with the above argumentation based on their 
angular spectrum Q , fulfilled by most optical wavefields 
outside near-held regions. In particular, this also involves 
that any far-zone scattered or radiated held, which is 
given by the propagating part of its angular ^ectrum [^, 
l67|, is gauge-invariant. This agrees with [5|,l9j and implies 
that since both £ and B are gauge invariant, the angular 
momentum J of these helds admits a decomposition J = 
L-l-^ into an orbital L, [proportional to ef • (r x V)M-I- 
p~^B ■ (r X V)C] and a spin [proportional to ef x M-l- 
p,~^B X C, see above], angular momenta; both L and ^ 
also being gauge-invariant. 

Therefore, it makes sense to study the conservation 
laws for the spin and orbital angular momenta separately. 
However, as we shall discuss in one example, sometimes 


Where the dissipation on transmission of spin how by 
interaction of the helds with matter is given hy = 
2 'k{£ ■ K, — j ■ C) and = (27rc/e/j,)(/C x B + J x A), 
respectively . 

Eqs. (IH) and (l6|) are conservation equations for the 
helicity and its how, respectively, and are as fundamen¬ 
tal as those for the energy, linear and angular momenta. 
We shall concentrate in (jh)) , where Si is linked to the ab¬ 
sorption by the particle, which gives rise, as shown later, 
to the mechanical action that it experiences as a conse¬ 
quence of the conservation of the helicity how. 

The time-harmonicity of the helds and potentials con¬ 
verts the quantities of Eqs. 0-0 into 


A = - 





( 7 ) 


and the time-averged (denoted by < • >): 


.3iS=<M’>= 




( 8 ) 


^ =< ^ >= —-S(eE* X E 
Ank 


1 


B* X B) + 

(9) 

and being the electric and magetic parts of ^ 
which now coincides with the spin angular momentum 
density. 

And 


-< >- 


CS[E*Bj + E*B, - ■ B],(10) 


which is the density of spin angular momentum how 
tensor. i,j = 1,2,3, 3 denotes imaginary part and 
k = nuifc; and we have written A = 5R[A(r) exp(—zwt)], 
C — 5R[C(r) exp(—iwt)], J = 5R[J(r) exp(—iwt)], 1C = 
5i[iir(r) exp(—zwt)]. Also in 0 and in 0 are re¬ 
placed by the time-averaged quantities: 


< 7r[^W^V • Cs{K x B*) 


ck 


-^3(J • B*) -f /^3(J . K* 
kn ck y e 


-K(E-K'*)]. (11) 
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< Mi >=-{cS(E X J*)i + 5R(B X K*), + 

■ K*) - ^{E,K* + E.K*)]}. (12) 

At this point it is pertinent to comment that in the case 
dealt here of time-harmonic fields, Maxwell’s equations 
and the above relations show that m and are pro- 
portional to Lipkin’s zilches [llj | , used in recent works as 
the field chirality and flow of chirality M , as 

well as the tensor of spin flow , 0: 

Mr=<je'>= = P <Mf> (13) 


y =<M >= = k^ > (14) 


Mij =< Mij >= k^JKj = k^ < jVij > . (15) 

The dissipative terms that appear in the continuity equa¬ 
tions for M and Mg, are however different from those 
of dU and dH). Ref. 0 argues that M and Jiij are the 
quantities with dimensions of angular momentum, rather 
than the zilches; and < >, < M > and < jVtj > are 

those to deal with in mesurements. We follow this, al¬ 
though for these monochromatic fields both kind of mag¬ 
nitudes should, as shown in (fl^ - (p® ■ lead to the same 
interpretation of effects. 


SCATTERING AND THE CONSERVATION OE 
SPIN ANGULAR MOMENTUM 

We integrate m in any volume v that contains the 
object with its distribution of currents. The divergence 
term is transformed into the flow of an entity across the 
surface E of this volume that will eventually be taken as 
a large sphere. 

Let us first consider a monochromatic, elliptically po¬ 
larized, plane wave incident on a body of arbitrary form 
enclosed in v. (cf. Fig. 1). The field at any point of 
the surrounding medium may be represented as the sum 
of the incident and scattered fields, thus its space de¬ 
pendent vectors are: E(r) = El®l(r) -|-El'*l(r), B(r) = 
BW(r) -pBl'’)(r). 

The incident field is written as e1®1 = 

Bb) = and that scattered in the far zone: 

El'*) = e(s) exp(ifcr)/r, BI'*) = b(s) exp(ifcr)/r . So that 
of course: b* = nSi x e®, e* • = b* • = 0; b = ns x e, 

e • s = b • s = 0. 

Then the time-averaged density of spin angular mo¬ 
mentum flow may be written: MCj = -p -p 



FIG. 1: (Color online). Illustration of an elliptically polar¬ 
ized plane wave incident on a polarizable particle of arbitrary 
shape. The fields are evaluated at the point P, of coordinates 
(r, 0, (j)) and position vector r = rs, (s being the unit vector 
along r), which eventually belongs to a sphere of integration 
of radius r, centered at some point ro of the particle, ro acts 
as the framework center 0. The point Q is the projection of 
P on the plane OXY ; the scattering plane being OPQ. We 
show the three orthonormal vectors: s, ey (in the plane OPQ 
and in the sense of rotation of 9), and e± (normal to OPQ). 

. Where 



2kn^JL ^ * 


-b 

eK 



-^yE(') 

*. b)* 



(lOa) 


2knfi * 


-b 

b*’ 




* . b)' 



(106) 



)* ^(s) 

-b 


'*Bf 


+Ef 


-b 

b'* 



)*. B^'*) -b E^'*) 

* -B^' 



(10c) 


From Eq.® we write the rate C at which this density 
flows on interaction of the incident wave with the body, in 
terms of its outward flux: £ = —{n? j4:TTc^) -s 

through the surface E of any large sphere of radius r with 
center at some point Tq of the body. dU is the element 
of solid angle and s denotes the outward normal. I.e., 
according to Eq.([6]): 

£ = £«-b/:'. (16) 

£(*), £(*) and £' are respectively the E integrals of the 
projections on s of — 

face of the sphere, respectively. Of course from Eq.®: 
£ = (n^/47rc^) drdflr^ < M >. 
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From these equations = 0, so that (HU) becomes 

£=£"+£'. (17) 

The integrals of — ■ s and —■ s across E leave the 

following non-zero terms, [cf. (fTOl) ]: 


£(«) = _ 


47rc^ 




SnkiJ 


Snkii 

1 

Sirkfj, , 


Js 

j dm{s ■ [e(">(s) X e(")(s)]}s = 

(in3{s • [b*''*^*(s) X b*'®^(s)]}s = 

/ dn3{es • [e^'*^*(s) x e*'®^(s)] 
JY. 


IbTrfc 


+ is. [b(")*(s) X b(*)(s)]}s. (18) 


and 


£' = 


47rc^ 


dVtr^^dsj = 


Snkfi , 


-f E(®)* • BW]s = 
K[B«*(ro)-e(s,)]s, = 


2^k? 


2^fc2 


5J[E«*(ro)-b(sO]s.. (19) 


In deriving (11911 we have used Jones’ lemma based on the 
principle of the stationary phase 6^, 7fll |: 




27rz 


[F{s,} 


— ikr 


-F(-s,)e*n- (20) 


Notice from and (ITU)) that all these quantities of the 
conservation law (nzi) - (HU) are also time-averages. Eq. 
(ITUl) conveys the ellipticity of the electric or the magnetic 
field, as well as the dual-symmetric one containing both 
e and b. This electric-magnetic duality is also seen in 

(HD). 

Ea. dTTll . together with (ITSll and (ITUll . constitutes the 
law for the moment of the force exerted by the incident 
wave on the object, expressed in terms of the flux £ of 
the spin flow tensor across any closed surface sur¬ 
rounding it. As seen, this stems from the conservation 
of the spin angular momentum density and, like the op¬ 
tical force [s^ it has a contribution £('*^ of the scattered 
field plus one C arising from the interference between 
this scattered field and the incident one. 

This force moment £ constitutes the contribution to 
the electromagnetic torque on the scatterering object 
(acting the ellipticity of the scattered field as a factor 


in the integration) by extinction of the incident spin an¬ 
gular momentum, and its transfer to the body at a rate 
proportional to a certain amplitude component of the 
scattered wave, projected along the polarization of the 
incident wave as a result of their mutual interference in 
the forward direction s = s^. 


CONSEQUENCES FOR A 
MAGNETODIELECTRIC BI-ISOTROPIC 
DIPOLAR PARTICLE 

Incident plane wave 


Let us consider a magnetodielectric bi-isotropic parti¬ 
cle, dipolar in the wide sense, i.e. characterized by its 
polarizabilities: electric , magnetic am, and magneto¬ 
electric Uem, ctme, such that for example if it is a sphere 
these quantities are given by the first order Mie coef¬ 
ficients oi, bi and ci as Oe = i-^cii, otm =_LSzki, 


ctem = ame = i^di, respectively [M|, 
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Then the electric and magnetic dipole moments, p and 
m, on the particle, induced by the incident field are: 

p = aeE«+aemB«, m = a™eE«+a„BW. (21) 

If, in addition, such sphere is chiral, then ci = —di 
and hence Oem = —ctme- 

The amplitudes of the fields scattered by this particle 
in the far-zone are: 

e(s) = k^- —[e“^(s x p) x s — ^ m)], (22) 


b(s) = k^- —[/i(s X m) X s -I- ■W^(® P)]- 


(23) 


Introducing (|22]) and (l23l) into (|T8)) and (ITUl) . evaluating 
the angular integrals, and substituting the results in (HZl, 
we obtain: 

k\ 


£(®) = ——5{-p* X p -b ^m* X m}. (24) 

6 e 


and 


£' = -^5R[(ip.B«*- 
e 


fim 


•E«*)s,]. (25) 


Eg. (1241) is half the rate of radiated angular momentum 
by the dipolar object [68| |. or equivalently half the re¬ 
coil optical torque on the particle which one would ob¬ 
tain from the conservation of the total angular angular 
momentum by inte grat ion of Maxwell’s stress tensor : 
J^dRdilrs X TijSj 371. 

Interestingly, since p • = m • = 0, (1251) may also 

be expressed as: 


£' =< To > . 


(26) 
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< Fq > is the extinction electromagnetic torque^ anal¬ 
ogous to the extinction energy in the optical theorem 

'kum, 


< Fo >= -5i[p X 


m X 




(27) 


Sometimes < Fq > is called ’’the intrinsic torque” |32l . 
47| exerted by the incident plane wave on the particle. 
However, Ea. d^ . is seen to be also an intrinsic 

torque component. 

It is useful to consider a Cartesian framework (cf. 
Fig.l) where the elliptically polarized incident plane wave 
has Si = (0,0,1). Its electric vector being written in 
an helicity basis = (l/-\/2)(l, ±z, 0) as the sum of 
a left-hand and a right-hand circularly polarized plane 
wave, so that Si = {eix,eiy, 0 ) = e'le^ + 


e, e 


and 

bi — (^ixj^zyjO) — ^( ^iyjCzTjO) — ^i ^i C — 

—ni{efe'^ — e~e~). The sign ± standing for LCP (-I-) 
and RCP (-), respectively. In such representation, the in¬ 
cident helicity acquires the form: = {e/k)^[e*^eiy\ = 

{e/ 2 k)S 3 = {e/ 2 k)[\e'l\'^ — |e“p], which clearly expresses 
this magnitude as the difference between left circular and 
right circular intensities of the field, extensively exploited 
in dichroism [l^, [13, llll and enantiomeric molecule dis¬ 
crimination. S 3 = 23[e*2.eiy] = |e)*~P — |e“p represents 


the fourth Stokes parameter [3, llC 

In the case of circular polarization of the incident wave, 
Cix = e, e is real, and Ciy = ±fe depending on whether 
the light is left or right circular; also nei = the 

incident helicity is = ±e^/fc. In addition, when the 
particle satisifies the first Kerker condition 71 , thus 

yielding zero differential scattering cross section in the 
backscattering direction, a^je = /rpm, then p = ±mm 
and b(s) = =Fme(s); i.e. the scattered field is circu¬ 
larly polarized with respect to the Cartesian system of or¬ 
thonormal vectors: (ex,e||,s), (see Fig.l), e± and e|| be¬ 
ing, respectively, perpendicular and parallel to the scat¬ 
tering plane |55l |. 

In consequence, in left circular polarization one has 
that e(s) = (e(s) • ej_)(l,±f,0) and b(s) = (ne(s) • 
e-L)(Tbl,0). The helicity of the scattered field being 
proportional to its intensity: = ±^|e(s)p, so that 

since then the flow of helicity becomes proportional to 
that of energy, one has that the torque £ , Ea. dTTi) . and 
the flow tensor ,Eq. OT, become proportional to the 
time-averaged force < F > and Maxwell’s stress tensor 
< >, [silTi: £ = ±(l//c) < F >, ±(1/^) < 

Tij >. This matches with a statement of [l3| on the ex¬ 
istence of a mapping of the helicity into the energy for 
circularly polarized waves. 


Arbitrary incident wave 

The above equations are generalized for arbitrary illu¬ 
minating wavefields, expressible in terms of its angular 


spectrum of plane wave components 66|, l67| : 

E«(r)= [ 

Jt) 

B«(r)= [ b«(s)e*'=(®-'’)dll. 

Jv 


(28) 


The integration being done in the unit sphere £> whose 
element of solid angle is dVt. In this regard, notice that 
Eas. (l22)l and (l23ll now constitute the angular spectra of 
the electric and magnetic scattered vectors when they 
are represented in the form (1281) . This integration con¬ 
tains both propagating and evanescent waves [illiilill; 
and to include them both, in (|2^ must be replaced 
by s*, complex conjugated of Si = (sf,s^,sf), where 

sf = ^1 — (sf'^ + s)'^) if < 1, (propagating com¬ 
ponents); and sf = i^J (sf^ -|- sf^) — 1 if sf^ + > 1, 

(evanescent components) . 

Erom (|28l) , using the same procedure as before for each 
plane wave component, and summing up all of them, an 
easy calculation shows that remains as in (1^^ . how¬ 
ever the extinction torque £' now becomes 


n 

2 kfi 


£' = - pmjdiEf’*}] 


?{i)* 


(29) 


which may be also expressed as 


1 


.1 


Where < Fq > now is 


-//(m • V)eW*}(30) 


< Fo >= isRip X E(*)* 


m X 




(31) 


The fields E^*) and being given by Eqs. 

At this stage it is worth observing how taking (I^Tl) into 
account, Ea. d^ confers, after separating real and imagi¬ 
nary parts, a physical meaning in terms of the extinction 
torque £' to some quantities associated to ’magnetoelec¬ 
tric effects’ introduced in [76|. 


On the torque involved in the conservation of the 
spin angular momentum 

The above equations not only show that the £ quanti¬ 
ties have torque dimensions, as they should, but as men¬ 
tioned before they manifest that £® is also an intrinsic 
torque. Notice on comparing (15(11) and (ISTl) with (1551) and 
(l27l) . the two additional terms in (l30l) due to the spatial 
structure of the incident field when this is not a plane 
wave. This term accounts for an extrinsic torque compo¬ 
nent. Also, the second and third terms of (|30p have an 
analogy with those of the dipolar component of the elec¬ 
tromagnetic force, whereas those contained in < Fq > 
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keep it with those of the Lorentz’s component of the op¬ 
tical force. This duality between the E and B vectors 
is also evident by comparing (Ei) with the expression 
of the electric (e) plus magnetic (m) forces on a mag¬ 
netodielectric dipolar particle In this connection, 

the recoil optical torque, Ea. dMl) . due to interference of 
the scattered helds and which, as said, is also valid for 
an arbitrary incident field, has a formal analogy with 
the electric- mag netic dipole interaction electromagnetic 
force: 01- 


In summary, the torque on the particle £ = £*-!-£', 
Eq. (IT71) , governed by the conservation of the spin angular 
momentum, is given by the sum of (l24l) plus (l2^ or (l30t . 
This torque has the recoil component added to £', 
which is as necessary to describe the dynamics as the 
above referenced Fe_m component of the force . 


£' -|-£® , Eg. ITTll . First, £', Eqs. (l29t or (l30)) , becomes: 


X 5(E X B*)] 

+ - ^Vx < S > +V5?(E • B*)] 

+ < S > -V»(E . B-)| 

-87r-(a^gV <We> < Wm >) 

€ 

2nk // ^ ^ 


t 

C V C ^em) ^ ® 


However this scattering torque (IMll is only half the re¬ 
coil torque exerted by the field on the particle, obtained 
from the conservation of total angular momentum 371. 
As we shall see, the other half is obtained from the con¬ 
servation of the orbital angular momentum. As expected 
from the premises above adopted in this formuation, this 
is consistent with the addition of both angular momenta 
to describe the total transfer to the body. 


DECOMPOSITION OF THE TORQUE 
INVOLVED IN THE CONSERVATION OF SPIN 
ANGULAR MOMENTUM INTO CONSERVATIVE 
AND NON-CONSERVATIVE COMPONENTS 


Like the electromagnetic force M, ^ M\, the optical 
torque is amenable of a decomposition into conservative 
and non-conservative components, as shown next. 


At this point we emphasize again that the vectors E*^®) 
and B*^®) satisfy the transversality condition V • E^®) = 0 
and V • B*^‘^ = 0. From now on we omit the (?) su¬ 
perindex. Taking in (1211) the real and imaginary parts of 
the polarizabilities, and using the vector identities: 


5R[(B* • V)E] = 1{Y + X); 5R[(E • V)B*] = ^{Y - X). 

3[(B* • V)E] = + R); 3[(E • V)B*] = ^(5 - R). 

A = 5i[V X (E X B*)]; 

Y = 5R[V(E • B*)] - kn^[F* x E] -h ko^[B* x B]; 
i? = S[V X (E X B*)]; S = 3[V(B* • E)]; ko = uj/c, (32) 


we obtain the following expressions for the parts of £ = 


The helicity Jif and the spin angular momentum density 
that appear in (13311 where defined in Eqs.ljH]) and 
respectively. Moreover, < We >= (e/167r)|Ep and < 
Wm >= (l/167r/i)|Bp are the time-averaged electric and 
magnetic energy densities, and < S >= (c/87r/r)3?(E x 
B*) represents the time-averaged Poynting vector. The 
R and / superindices of the polarizabilities denote their 
real and imaginary parts, respectively. 

Notice that there are conservative terms containing 
af^koXJff and a^koVJ^f, which represent a gradient op¬ 
tical torque where now the helicity JY' plays a role anal¬ 
ogous to that of the field energy for the optical force. 
This may constitute the basis of what we may call as 
a rotational optical tweezer, positioning the particle in 
equilibrium points where there is no lateral torque. This 
is addressed in some of the examples below. Also the 
imaginary parts and yield gradient terms of the 
energies, like in the gradient force components. These 
latter terms of the cross-polarizabilities change sign if the 
chirality of the particle, characterized by aem and Ume, 
is reversed, a fact that suggests its use for enantiomeric 
separation [^ . Likewise, there are conservative terms 
V5?(E-B*), this time with the imaginary parts a], and 
■ 

The polarizabilities and also appear in curl non- 
coservative parts that contain the vortex circulation of 
the reactive alternating flow of energy V x 5(E x B*) [ei]. 
These latter terms are sometimes zero, like for a highly 
collimated circularly polarized Gaussian beam. On the 
other hand, the a^, and non-conservative terms pro¬ 
portional to come from < Eq > in Eq. dsni). They are 
the analogous to the radiation pressure terms containing 
the held linear momentum in the time-averaged force. It 
is remarkable the analogy according to which while in this 
latter case the conservation of the wave linear momentum 
leads to the radiation pressure part of the optical force, 
in the optical torque the conservation of the spin angular 
momentum ^ leads to terms that similarly contain this 
quantity and play the role of a spinning torque. 












There are also energy flows in the cross-polarizabilities 
a^g and afJn, that contain the curl of both the Poynt- 
ing vector and the flow of spin angular momentum, and 
thus represent their respective vortices surrounding the 
particle. 

Likewise, the recoil torque, Ea. d^ . is written as 


„ f 4:Tlk r 1 I 19 9 I I 9 

C = -—{-am 

6 c 

+ “T|aemk^m + 

- fj-^{a^ame)] < S > 

c e 

-2[^^{a:aem) + x B*)}. (34) 


which in addition to non-conservative terms with elec¬ 
tric and magnetic spins, and with the energy flux, has 
conservative terms, if the particle is bi-isotropic, contain¬ 
ing the reactive or stored energy and its alternating flow 
[68l| : 3(E X B*), which taking into account the identity: 
A:3{E X B*} = (1/2)V|E|2 - SR{(E* • V)E, yields gradi¬ 
ent force-like terms, proportional to V|Ep, which change 
sign if the particle chirality varies from left-handed to 
right-handed. 

Thus in summary we see that the helicity and the spin 
play a role in the optical torque somewhat analogous to 
the energy and the Poynting vector in the optical force. 

In the particular case in which the particle is chiral, 
then [l^ aem = —ctme and (I55)) - (I5S1) reduce to 


X 5(E x B*)] 

+ - ^Vx < S > +V5?(E • B*)] 

2e c c 

X < S > -V5?(E • B*)] 

2 c c 

„ M r „ 2nfc u p „ _ 

-STT-a^gV <w> -X ^ 

e c e 


IGirkfi fjl / c 

r^me < S >}■ (35) 


o k ^ 4:Tlk , 1 I 19 _ 9 I 19 

C- = -—{ ^ [^jOel =9^6+1* IcTml •^m 

-k\0ime\^ + M^'^e)] 

-^^^[-3(Q;*ame) + n^ia^ame)] < S > 
c e 

+2[^?ft{a:a^e) - M5R«ame)]3(E x B*)}. (36) 


become 


n ^ 2ai ,k^ ^ 27r^ „ , 

£' = *oV^H-^ ^-Vx < S > 

k € 6 C C 


1 j 

—ATT-aL,^V < w > -- 


nk 1 


X ^ 


Snk ni ft „ T ___ 

- —< S >}., (37) 


= 


2fc% 


-87r^3(a:a^g) < S >}. (38) 
e 


Thus the reactive parts dissapear. We obtain two kind 
of terms in (1571) and (155)) : those conservative in C linked 
to the gradient of the helicity and energy, and those non 
conservative in C and £® containing the flows of helicity 
and energy, as well as their curls. Once again, this has a 
suggestive analogy with the gradient, radiation pressure 
and curl components of the optical force 541. 


THE ORBITAL ANGULAR MOMENTUM 

In the non-absorbing dielectric medium of refractive 
index n = y/ejl the field density of orbital angular mo¬ 
mentum L and of its flow Ay, (i, j = I, 2, 3), are defined 
in the dual-symmetric formulation here employed, as Q : 

'L = ^[£-^A + B-^C], (39) 

if = r X V. (40) 


+BjAi-£jCi], (f,j,fc,Z, TO, n = 1,2,3). (41) 


eiki is the antisymmetric Levi-Civita tensor. In (1391) and 
(1511) A and C are the vector potentials introduced in Sec¬ 
tion 2. 

Under the premises established in Section 2, L and Ay 
fulfill the continuity equation Q 

L -k V • Ay = (42) 


Where represents the dissipation by transfer of or¬ 
bital angular momentum on interaction of the fields with 
matter. 

In addition to o, the time-harmonicity of the fields 
leads to the time-averages Q: 

L =< L >=-^S{eEU^E-k ^B*-.ifB}, (43) 

4nK /r 


< w > denotes the density of electromagnetic energy: 

< w >=< We > + < Wm >■ 

It is interesting that if the particle, in addition to being 

tte = then dSlD and OMD 


chiral, is dual 71j 


mg tnat 

zS-Izi: 


^-1 


Ay -< Ay >- 




+B*E,+B*E,}, (*,j = 1,2,3). (44) 
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Where in (l42l) is replaced by its average: 

'KC 

<J(>= — {3[J* • ^E] - • ^E] 

nk 

+ E*K,]}. (45) 

Xk being the fct/i-Cartesian component of the position 
vector r. 


CONSERVATION OF THE ORBITAL ANGULAR 
MOMENTUM IN SCATTERING 


Equation (l42l) expresses the conservation of orbital an¬ 
gular momentum flow and, again, consveys a density of 
force moment, or torque. We next find what does this 
mean. As before, we integrate it on any volume sur¬ 
rounding the scattering object. The divergence term is 
then transformed into a flow across the volume surface E 
that will eventually be taken as a sphere of large radius 
r. 

With reference to Fig.f let us first consider a 
monochromatic, elliptically polarized, plane wave inci¬ 
dent on an arbitrary body. The field at any point of the 
surrounding medium is again represented as the sum of 
the incident and the scattered field. Then the flow (or 
time-averaged flow) density of orbital angular momen¬ 
tum may be written by means of A^ decomposed as: 
Ay = Ag^ + A^g + A'j- . Where 


i{i) 

Hj 


c2 


*Ef + Bf *Ef'^] 


')}. (46) 

(47) 






-Eg 




+bI"'> *^E^^ - 


From Eq. (1421) the rate at which the total or¬ 
bital flow varies on interaction of the wavefield with 
the body, is given as before by the outward flux: 
— (n^/47rc^) f^dUr'^Aij ■ s through the surface of any 
large sphere S of radius r with center at some point Tq 
of the region occupied by the object. I.e., according to 

Eg.dia: 


= (49) 

Where 4)^, 4)^^ and 4)^ are the S-integrals of the pro¬ 
jections on s of — A-g, — A-®^ and — A'^, respectively. 
4*.Sf = (n^/47rc^) J^drdEtr^ < ^ >. 


From these equations we have that 4>y = 0, so that 
dUl) becomes 


= 4>^) 


4>', 




(50) 


Operating with (|i5|) - (|i5|l and using the identity: 
^jik^kmn — dj^Sin the integrals of A 


(«) 


and —A',- • s across S lead to: 


47rc^ 

1 




— —- / d0r^3[ee*(s) • Jjfe(s) = 
oirk Jy: 

b*(s) • ^b(s) = 


Stt/c 


f d0r^5[ee*(s) • .ife(s) -I- 
IOttk Jy 


1 


b*(s) • Jfb(s)]. 


(51) 


And using again Jones’ lemma, Eq. (EHl), we get: 
^2 r 


1 


dOuAySj = 


__5i{eeW*.[^e(s)],=,J = 


(52) 


Eqs. (ISTl) and (l52l) constitute the conservation law for 
the orbital angular momentum of either the electric or 
magnetic field, as well as of the dual-symmetric one with 
both e and b. Fa. (150)1 . together with (1511) and (I52p de¬ 
fine the force moment, or torque, by transfer of this or¬ 
bital angular momentum from the incident field to the 
object, expressed in terms of the flow of the tensor Ay 
across any closed surface surrounding this body. It is 
very interesting that the rate at which this incident or¬ 
bital angular momentum Eg. (1501) . is extinguished 
by transference to the body is proportional to a sort of 
cross-orbital angular momentum given by the interaction 
of the incident field with the scattered field in the forward 
direction s = Si. 

Moreover, recalling that the incident field is a plane 
wave, and since r = rs, one has for the incident electric 
orbital angular momentum 

cj[E(d* . = ^{E^^>ikr^[s, ■ sjE^*)} = 

S{e^*^*iA:r[s x = |e*'*g^A:r(s x s^). (53) 

With an analogous expression for the magnetic momen¬ 
tum. Then the total incident orbital angular momentum 
enclosed in a sphere of large radius i? is s given by 

pR p'k 

/ drr’^ / d(l) desin0^[F,^^'>* ■ 

Jo Jo Jo 

= 0 . 
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As expected. Therefore the extinction torque from this 
orbital angular momentum, Eq. (I52L is zero if the in¬ 
cident wave is plane. (However, as we shall see next, in 
general it is not zero for any arbitrary field). 

On the other hand, the physical meaning of the flow 
of Eq.®, corresponds to a force moment of the 

scattered field on the particle. This is the consequence for 
this recoil torque of the conservation of orbital angular 
momentum. But since since the incident plane wave has 
no orbital moentum, this recoil contribution has to come 
from a transfer of the incident spin angular momentum 
into the scattered orbital angular momentum. Hence this 
effect constitutes the essence of the spin-orbit interaction 
Si as a consequence of the scattering of the incident 
field by the object. 


EFFECTS ON A MAGNETODIELECTRIC 
BI-ISOTROPIC DIPOLAR PARTICLE: THE 
ELECTROMAGNETIC TORQUE INVOLVED IN 
THE CONSERVATION OF THE ORBITAL 
ANGULAR MOMENTUM 


Incident plane wave 


Let us now address a magnetodielectic bi-isotropic 
dipolar particle. Using Eqs. ®-(l23l), evaluating the an¬ 
gular integrals of eiD, and substituting the results in 
dSQl), we obtain the extinction torque from the orbital 
angular momentum transfer to this scattering body: 

(54) 


Likewise, we get the scattering, or recoil, torque from the 
orbital angular momentum: 


^3 1 

= —^3{-p* X p -I- /rm* X m}. 


(55) 


Eq. (1551) is identical to (l24l) for the scattering torque from 
the spin. Both recoil torques sum to yield the resul¬ 
tant scattering torque, which coincides with that derived 
from the conservation of the total angular momentum J 
through the integration: J^rs x TijSjds, [37|. 

On the other hand, Eqs. m and the transversality 
of the incident field show that given by (l54ll is zero, 
as expected from the above discussion. However Ea. (|55ll 
exhibits the non-zero scattering contribution (recoil 
torque) to <i>^. Eg. (1501) . which as mentioned before ex¬ 
presses the spin-orbit interaction involving the conversion 
and transfer of the spin into orbital angular momentum 
torque by scattering with the object. 


Generalization to an arbitrary incident field 

The above equations are generalized to an arbitrary il¬ 
luminating wavefield expressible by its angular spectrum 


(l28t . which introduced in (l54t yields 

4>^ = ^y|s{i(p-V)B*-Mm-V)E*}. (56) 

This form is equal to half the contribution from the spa¬ 
tial structure of the incident field to the extinction torque 
derived above from the spin conservation, as seen by com¬ 
paring (1551) with the second term of (1501) . 

On te other hand, Eq. (1551) for remains valid for 
an arbitrary incident field. Thus we obtain from (|50l) 
using ([551) and (I56p the following torque involved in the 
conservation of orbital angular momentum: 

4/t Y ^ € 

1 

-3{-p* X p -I- ^m* X m}. (57) 

6 e 


DEGOMPOSITION OF THE ANGULAR 
MOMENTUM TORQUE INTO CONSERVATIVE 
AND NON-CONSERVATIVE COMPONENTS 


Using (I0T1) and taking into account the vector identi¬ 
ties (1551) . we express the two components (1551) . and 
d)^, (1551) . of the orbital conservation torque, (1571) . of the 
transversal fields, again droping the superindex (f), as 


+a^[kop.‘^VJ^ + X 5(E X B*)] 

+ - ^e) - —Vx < S > +VK(E • B*)] 

2c c c 

_ ^^) _ ^vx < S > -V5?(E • B*)] 
2 c c 

-87r^(a^gV <We > < Wm >) 

2nk jjL 


t 

+ ^(aL + aLmExB*)},(58) 


and 


„ k^ Ank , 1 , ,o _ _ 

4’^ = —{ ^ [^loel Ae + fj- \oim\ Am 

-\--^\otem^'^m + | Ome 

+ ^^p^[^S(a*aem) - p,^{a*^ame)] < S > 

-2[^ift{alaem) + (59) 







11 


If the particle is chiral, aem = —ctn 
duce to 


= ^{a^[ko-Vjr - X 3(E X B*)] 

rCIJj 6 ^6 


and (l5^ re- 


Where the electromagnetic extinction torque is 

< r' >= £' + = 

4A: Y /i e 


+a^[A:o/x^V^+ x ^(E x B*)] 

+ - ^e) - —VX < S > +VK(E • B*)] 

2e c c 

_ ?!I£vx < S > -V»(E. B-)] 

2 c c 

-STT^a^gV <w> -^^a^eV x ^}. (60) 


o k f ^nk r 1 I 19 91 19 

=-^ M jcTml 

-k\0ime\^ {-^^m + 

-^^^[-S(a*ame) + ^l^{a*^ame)] < S > 

C € 

+2[i5R(a:a,^e) - M5R«an9e)]S(E x B*)}. (61) 

And when the particle is dual, (ae/e = /iCTm), Eqs. (ISOl) 
and (El) become 

n cv^ 

$' = {^fcgVjr 
fe e 

2ai ,k‘^, ^ ^ , 

H-[- {i^m ~ -^e) -Vx < S >] 

e c c 

-47r-aLV < u; >- —-a^V X ^}. (62) 

e c € 




And the electromagnetic recoil or scattering torque reads: 

< r" >= c 


k^ 1 

= ——3{-p* X p + ^m* X m}.(66) 
o c 


0^3 1 

$S, = -^{nfc[^|aep^ 

+ \0lTne\'^ { —+ M^'^e)] 

-S-K-^{alame) < S >}. 


(63) 


Similar remarks as for the analogous decomposition of 
the torque from the spin, above, apply to these latter 
equations. 


THE TOTAL ELECTROMAGNETIC TORQUE 

From all the previous analysis we conclude that the to¬ 
tal electromagnetic torque on the particle, < F >, stem¬ 
ming from the conservation laws of both the spin and or¬ 
bital angular momenta, is given by the sum of £, Ea. (ll7L 
[cf. also (l24l) and (l29l) or (l30l) ] and Ea. dS^ . which 
we express as: 


< r >=< r' > -I- < r® > . 


(64) 


As mentioned above, we emphasize that this total torque 
is also obtained from the conservation of the total angu¬ 
lar momentum < J > 37|. The analysis given here, how¬ 
ever, elucidates the contribution of both the spin and the 
orbital angular momenta through their respective conser¬ 
vation laws, and shows that both torques are added, like 
the angular momenta. 

We notice that making use of El one has, (dropping 
the superindex (i) in all following equations, understand¬ 
ing that we deal with incident fields), 


< Fo > -f < F'* >= 


Qfg fc^Jofel 

T ~ 


+k-\ 


,al, k^ , , 


\a 


a„er)]3(E* X E) 
2 

-)]5(B* X B) 


+ [ 


_ O/^ 1 

+M3(a^a;,J]5R(E x B*) 
-4 1 2fc3 1 


me ' em 


-—(-??(«: aer.) 

z, o e 

+MK(a^a;,J]5(E x B*). 


(67) 


However, introducing El in the optical theorem that 
represents the conservation of energy M | 


rk‘^ h! 

= 2^[P• < + • K]; (68) 

W“ standing for the rate of energy absorption by the 
particle, one has 


LU 


,a. 


a{ fc^Jdel 

"2 3"^~r 

fc3 


2- + 


+ lJ^\0ir, 


+[ 


Oi^ — Oi^ 
^me ^em 


2k^ A 


i]|E|^ 

—;]|Bp 

^{alaem) 


+/iS(a^a;,J]3(E • B*) 


^me ' ^em 


2fc3 1 


R(a*aem) 


+^iifiia^a*^,)]iR(E ■ B*). 


(69) 


It is interesting that for fields such that B = =fE, (like 
e.g. the circularly polarized field scattered from a dual 
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particle, mentioned in Section 4), one has from (I67|) and 
(l69)) that 


< To > + < >= X E) = ^^e.(70) 

Where we have defined the absorption cross-section as: 

= (87rA:/a;|Ep)>V^“^. This occurs in particular for 
an incident circularly polarized plane wave and for some 
beams, as well as for a pure electric, or a pure mag¬ 
netic, dipolar particle, [in this latter case #e should be 
replaced by in (1701) ] . In such situations we see in (1701) 
that the component < Eq > of the extinction torque 
< r' >.[cf. (10^ ]. is cancelled by < E® >, thus only re¬ 
maining in the sum (1701) that part of < E® > which con¬ 
tains . 

More generally, for any wavefield, if the particle is mag¬ 
netodielectric, although not bi-isotropic, (aem = Ome = 
0), separating in (|69l) the electric and magnetic parts by 
writing W“ = W“ -|- and introducing the electric and 
magnetic cross-sections: = (87rfc/a;|Ep)VVi“^ and 

= (Sirk /one has on writing the optical 
theorem (1001) separately for the electric and the magnetic 
parts. 


< Eo > -f < E® >= x E) 

Sttk 

+ ^ ^ (71) 


These equations are seen to be compatible with 
(l36t . (l58l) and (l59l) . 

Hence, in contrast with some previous work 27-^ 
and in a gree ment with some experiments and calcula¬ 
tions [lil [Sfl [ssj, these equations show that 

this particle does not experiences a torque and spins due 
to the so far called intrinsic torque: < Eg >, [which, 
as we have just seen, is cancelled by a part of the re¬ 
coil torque on this kind of particles through the optical 
theorem relating imaginary parts and moduli of the per¬ 
mittivities, Eg. (15^ ]. but because it absorbes the incident 
energy and, as shown by dZOD and (ED, receives the spin 
or through the remaining part of the recoil torque 
that contains the absorption cross section. On the other 
hand, as we see below, the orbiting of the particle is due 
to the beam shape carried out in the torque by the second 
term of (p51) . 

In this regard, it is thus pertinent to remark that in 
contrast with previous torque theories that do not com¬ 
ply with the optical theorem, the proportionality factors 
between the optical torque and the electric and magnetic 
spin angular momenta is given by the absorption cross 
sections, not by the imaginary part of the corresponding 
electric and magnetic polarizabilities. This is a conse¬ 
quence of the static starting point and hence of the lack 
of consideration of the recoil torque in such theories. 

It is well known that the static regime applies for 
Rayleigh particles, (for which ka << I, a being their 


linear dimensioUjLe. for example their radius if they 
are spheres) (m. 57. US- In this limit the optical torque 
from such theories coincides with the one obtained in this 
paper. 

If for example the dipolar particle is purely dielec¬ 
tric, (am = ctem = OLme = 0, is just denoted as 
cr(“)), Eq. dM)) which becomes —(2/3e)fc^|ae|^ = —a], + 
£cr(“)/47rfc, introduced in the recoil torque < E® >= 
(l/3e)|ae|^S[E* x E], [cf. Eg. (1571) ]. leads to cancella¬ 
tion of < Eg > by the corresponding af term of < E® >. 
Hence, Eqs. (l64l) . (|65]) and (1671) lead to a torque on this 
particle: 

gcr(“) 


< E >= 


Stt/c 


■3[E* 


X E -k 


4fc 


"3{i(p-V)B*}. (72) 


We see at once that the first term of (1721) . which we ex¬ 
tend to right side of ED, and that represents the spin 
torque that we denote as < E^- >, characterizes the rota¬ 
tion of the particle on its axis due to the spin transfer by 
absorption of the incident energy through . 

In connection with the aforementioned static ap¬ 
proximation, notice again that the above relationship: 
— (2/3e)fc^|aeP = —af -k /Airk becomes: af = 


jAmk only in the Rayleigh limit [sj, 57, 0. 

Concerning the orbital movement due to the field 
structure, generally involved in the second term of (|65p . 
it is useful to consider the vast kind of wavefields that 
posses a vortex Expressing E(r) and B(r), as 

well as the position vector r, in cylindrical coordinates 
{R,(j),z) and extracting the vortex phase we write 
such incident fields as 


E(r) = e*'‘^E(r), B(r) = e*'‘^B(r). (73) 

I being the topological charge. So that from ED and 
after Eq. ED, we express this term of (1551) as 

^yi5{i(p-V)B*-/r(m.V)E*} 

= - I^«L)[S(E • V)B* - 

+ (ial - /.a«)[3(B • V)B* - ^K(H,^B*)] 

+ {-^ai - tiaUim ■ V)B* + 

+ - /ra^)[5?(B • V)B* + i^{B^B*)]}. (74) 

e K 

For example, if the particle is purely dielectric and 
isotropic, (only ae 0), then the right side of (1741) re¬ 
duces to 

Av/l9(l(p.V)B- 

+aimB-y)B* + ^^{E^B*)], (75) 
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which accounts for the particle orbital movement, and 
that expressing again al, in terms of through the 
aforementioned optical theorem: o;^ = /Airk + 

(2/3e)fc^|aeP, shows that there is an 1-number torque 
component, proportional to la^°‘'^/R, making the particle 
to orbit around the vortex axis via the transfer of orbital 
angular momentum by absorption; also existing another 
orbital component by transfer through the particle scat¬ 
tering cross section involved in l\ae\’^/R [5J, l72|. Both 
components are amenable of observation . 

For fields with a well defined incidence direction, like 
optical beams, the term (3/4fcn)ag5R(E • V)B* of (1751) . 
which plays for the torque a role analogous to the 
(l/2)ag3{(E • V)E*} component of the orbital momen¬ 
tum in the optical force 52, , may make negative the 
resultant longitudinal torque (i.e. opposite to the inci¬ 
dent helicity). This is illustrated below with incident 
Bessel and Gaussian beams. 

There are more components of this extrinsic torque, 
(ESI, in the transversal 4> and R directions; but for the 
purely dielectric particle the z component of the two 
terms of in (1751) is zero for a variety of beams, like 
those discussed in the following examples. More complex 
effects of this sort remain to be studied in Eg. (174)) when 


^em ^nd OLjne 


are non-zero. 


EXAMPLE 1: INCIDENT BESSEL BEAM 

As a first illustration we consider a transversal electric 
(TE) Bessel beam [i^, propagating along OZ in air, 
with a vortex of topological charge I, whose electric and 
magnetic vectors are 

E«(r) = ^e^^^e^^‘^{-^Ji{KR) +^^J^{KR), 

-^RiKR)-z^JliKR),Oy, 

BW(r) = 

- — jyKR)+i^j{{KR), j^MKR))-, ( 76 ) 

k — {kx , ky, hz'), K — : ky: ^ — |k| — ~, 

e = n = n = 1; r = {x,y,z),Ii= {x, y, 0), 

incident on a dielectric particle whose induced dipole mo¬ 
ment is p = agE^®). 

This beam is eliptically polarized. Eirst, we note that 
the helicity of this field is 

Ik K 

J^ = e^^Ji{KR)J[{KR) ( 77 ) 

Then according to (l29l) or (l30l) : 


With 


g2 Iff 

<To >=ai^ — Ji{KR)Ji{KR){0,0,l) 

= ^af^(0,0,l). (79) 

Which shows that this component of the torque follows 
the helicity. Also, Eg. (1751) manifests a contribution of the 
vortex charge I to this part of the torque. 

On the other hand. 


>C' = 2$^ = ^S{(p-V)bW*} = 

p^k IK , 

-^{afj^[KJi{KR)J['{KR) + KJ^^KR) 

-^Ji {KR)J[{KR)\{x, y, 0) + ai [[^ jf {KR) 

-2’^JiiKR)JliKR) + ^j;\KR)]{y,-x,0) 

IK^ 

-2 — Ji{KR)J[{KR){Q, 0,1)]}. (80) 


Where J[{KR) = dJi[KR)/d{KR). 

Eg. dSD)) exhibits an azimuthal dependence {y,—x) = 
—Ref) of the al terms, (i? and (f> are the cylindrical co¬ 
ordinates in the transversal plane of the beam). Hence, 
although < To > is longitudinal, i.e. along OZ, the beam 
structure creates and modulate a transversal azimuthal 
component of C. Moreover, there is a radial torque 
[x, y) = Rf in the terms due to the helicity gradient. 
This is extensive to the torque part: ^3{(p • V)B*^®)*} 
[cf. Eq. (1571) ] coming from the orbital angular momentum 
conservation law. 

It is interesting to analyze the different parts in 
Eqs. dZHl)-® from the point of view of the decompo¬ 
sitions (1551) and (1551) . There are the conservative compo¬ 
nents: 

1 Ih K 

-afV^ = —^^af[KjyKR)Jl'{KR) 
+KJ[\KR) - ^Ji{KR)J'i{KR)]{x,y,G) (81) 


^a^V5R(E.B*) = 0. (82) 

Eq. (ISTI) shows the aforementioned radial gradient torque 
whose sign oscillates with the distance R to the beam 
axis. 

The non-conservative components: 

-al^e = i < To > . 

c z 


( 83 ) 
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And 


Ic ^TT 

-ai^m - -^aiVx < S >= 




kc 

-2^^Ji{KR)J[{KR) + ^j[\KR)\{y,-x,Q) 
IK^ 

-2—MKR)JliKRM0,l)]. 


Having used: x‘^J"{x) + xJ'i{x) + (x^ — P)Ji{x) = 0. 
Also, 


- — afV X 3(E X B*) = 0. 


On the other hand, 


(85) 


k IK 

= -^e^\ae\'^ — Ji{KR)J[{KRmQ, 1); (86) 

which is contributed by the electric helicity flow term: 
of Eqs.(l36l) and (IMl) . 

Using as shown above the optical theorem: Og — 
(2/3e)fc^|Q;eP = j^-'Kk, and since + k1 = fc^, the 

resulting torque becomes 

p2 IT^ 

< r(r) >= —A(XA)J((i^i?)](0,0,l) 

QpZjL Iff 

+^{a^^[KMKR)Jl'{KR)+KJ^^{KR) 

-^MKR)JliKR)]{x,y,0)+ai{[^J?iKR) 

-2^^Ji{KR)J[{KR) + ^j[\KR)]{y,-x,Q) 

-2^J,(ifi?)J/(ifi?)(0,0,l)]}. (87) 

The spin torque < >, [see Ea. (l72]) ]. that results from 

the sum: < Fq > + 4>^, [cf. Eqs. (1551) and (1551) ]. 

is the first term of (EZD, and as said above, describes 
the particle spinning following the incident helicity by 
transfer of SAM through energy absorption. Eq. (1571) also 
shows that, in particular if the wave were plane, this term 
would be the only contribution to the torque, there be¬ 
ing no orbiting of the particle, as expected. Otherwise, 
there is a radial helicity-gradient term of a^, and an az¬ 
imuthal vortex term of in (1571) due to both the beam 
transversal structure and the interaction of the transver¬ 
sal E and B with the longitudinal B. The latter mani¬ 
fested through the particle scattering and absorption of 
the incident beam as seen from the optical theorem, and 
producing an orbiting of the particle around the vortex 
axis OZ at maximum intensity points, which are stable 
due to the gradient force. 

Eig. 2 shows three quantities in the transversal plane, 
where they have an annular spatial distribution, for a Si 
sphere, (refractive index np = 3.5) illuminated at a wave¬ 
length near its dipolar resonance [13,173, and if = 0.9fc. 



FIG. 2: (Color online). Doped Si sphere of radius 230 nm 
illuminated at A = 1350 nm, where its electric polarizabil¬ 
ity dominates, by a TE Bessel beam. K = 0.9fc, I = 2, 
of = 1.095 X 10''nm^ ai = 6.083 x 10®nm^ cr^) = 0.1. 
Red: Longitudinal component of the time-averaged optical 
torque < F >z. Blue: Longitudinal spin torque < Fc > 
transferred by absorption, (its magnitude appears multiplied 
by 10^). Broken pink: Helicity R of the incident beam, (its 
magnitude appears multiplied by 10®). The three quanti¬ 
ties have circular symmetry with annular spatial distribution 
in the transversal AF-plane, and thus in this representation 
they are even functions if the axis of the radial coordinate R 
is extended to the left of 0 . The torques and helicity are plot¬ 
ted in nm® and nm~®, respectively, since they are normalized 
to half the incident intensity density: e^. 


The particle is assumed to have been slightly doped to 
confer it a low absorption, so that its spin torque is ob¬ 
servable through The optical torque < F is 

dominated by the last term of (1571) . which is several or¬ 
ders of magnitude larger than the spin torque < F^- >. 
Also, while the latter follows the incident field helicity, 
the former opposes it, thus giving rise to a longitudinal 
z-negative optical torque on the particle, which conse¬ 
quently orbits around the vortex axis in opposite sense 
to its spinning. Diminishing K results in a weaker torque, 
but does not change its negative sign. Similarly occurs on 
varying /, which only changes the oscillation points along 
i?, slowly diminishing their ampliitudes as |^| grows. 

It is also interesting in connection with the gradient 
torque (1571) and the corresponding behavior of the helic¬ 
ity in Fig. 2, that this optical torque transversal com¬ 
ponent is maximum and zero at those points where Jif 
changes sign and is extreme, respectively. These posi¬ 
tions coincide with those where the gradient force mag¬ 
nitude on the particle is maximum and zero, respectively. 
Hence while this latter force tends to place the object at 
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an equilibrium position of minimum energy, the gradi¬ 
ent torque tends to a zero transversal value, both in its 
spinning and orbital movement, (cf. Eqs. (I55|) and (1551) 
for the spin and angular gradient optical torques). As 
anticipated in Section V, this effect can be exploited as 
an additional degree of freedom in what we may call a 
rotational optical tweezer. 

EXAMPLE 2: INCIDENT GAUSSIAN BEAM 


The ith component of the orbital momentum Pq being 

M 

hj = 1,2,3. 

With the mean energy flow: 

< S >= ^e2e“^(0,0,l). (93) 


TE and TM Gaussian beam 


And 


We shall now consider a wavefield commonly used in 
several works and illustrations. This is a Gaussian beam 
in air, both TE and TM, circularly polarized, with spot 
size 2cr and a vortex of topological charge 1: 


2 o2 O / 

S{(E* • V)E = (±^ - ^)(y, -x, 0). (94) 

Now observe that since for the paraxial TE and TM field 
(l88l) . V • E 0, one has 


EW(r) = e(l,±i,0)e*'V'=^e“^. 

BW(r) = 4-V X E«(r) ~ ^8, x E^^Vr) = z x E«(r) 
ik ik 

= TfEW(r) = e(Tb 1, (88) 

k — ( kx , ky^kz)^ K — ( kx , ky ,0), k — |k| — —, 
t = p = n = l\ r = (a;, y, z), R = (a:, y, 0), 


illuminating a dielectric particle whose induced dipole 
moment is p = OgE^*^. 

This is an interesting case because as we see below, 
although some aspects of the forces stemming from the 
above paraxial expressions of the E and B vectors are 
easily obtained, these fields are not divergenceless and, 
as we see next, are not appropriate for any theory based 
on the free-space Maxwell equations: V • E = V • B = 0, 
like the torque formulation developed here. Neither the 
total angular momentum density, nor its z-flux across a 
transversal AE-plane, can be expressed as the sum of 
the corresponding spin plus orbital densities and fluxes, 
as expected from the remarks of Section II. 

First, we note that the helicity of this field is 

1 rj p2 

M’ = ±—e^e~~^. (89) 

k 


It is convenient to analyze the dynamics exerted by this 
beam on a dielectric particle from the point of view of 


the mean optical force, which is |48l. l5ll. l52l. l54| 


<F +ai — Po. (90) 

4 c 

We omit the superindex (z), understanding that in all 
following equations we are dealing with the incident field. 
The gradient, pulling the particle to the beam axis, is 


4 2 

V|E|2 = —eVT^(-a:,-y,0). (91) 

{-x,-y,0) = -Rr. 


S{(E* • V)E = 3{(V • E)E*} - X S{E* x E}. (95) 

The existence of the non-zero first term in the left side of 
(1551) is a fact that should be recalled when expressing the 
orbital momentum in terms of the electric spin < > 

curl: V X S{E* x E} for fields that are not divergenceless. 

Notice that the orbital momentum Pq, described by 
(l92)) - (|94)) . contains the term ±^, [cf.Eq. (l94l) ]. which is 
essential to describe the spin curl component of the force, 
[cf.Eq. (1551) ]: 

1 4 2R^ 

-Vx3{E* xE} = Te"—e-T^(y,-x,0). (96) 

2 (j^ 

{y,-x,0) = -R^. 


Which together with 

S{(V • E)E*} = e^e-^ (t^ - ^){y, -x, 0), (97) 

yield the expression ((94l) when they both are introduced 
in (1951) . In fact, the omission of the leads to the 
well-known fundamental paradox of having a zero spin 
momentum, [the last term of Eq. (1951) 1 while the spin 
angular momentum is not zero [79|. 

Eas. (|90)) . and (l92l) - (l94l) . characterize the force due to 
Po, pushing the particle along z with the radiation pres¬ 
sure characterized by < S > , and in the azimuthal direc¬ 
tion due to 3{ (E* • V)E given by the divergence and spin 
curl terms of (1951) . Also, from the above equations one 
has that the total angular momentum J does not hold 
J = L -f where J = rx < S > and L = r x Pq. 
Neither their corresponding flows across the transversal 
AE-plane fulfill such decomposition. 

One would write the torque R x Og < S >, given 
by the moment of the radiation pressure with respect to 

the beam axis, as azimuthal: — e a^cp. Whereas 
the torque R x f^S{(E* • V)E from the moment of the 
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remaining orbital part of the force would be longitudi- 
2 

nal: —— ■^)z, (notice this orbital z- 
component being proportional to the vortex index I and 
independent of R), to which the moment of the spin con- 

j 2 

tribution R x ^(V x 3{E* x E}) is ±e^^e“"^ a^R^z. 
However, as we have already seen, the torque by transfer 
of spin to the particle actually comes from its absorption 
cross section. If there is no absorption, (or birefringence 
(^1. there is no spinning of the particle due to the SAM. 
On the other hand, according to (1^ or O: 

+ =< To > +2$(s.. (98) 

A V fj/ C 

With e = ^ = n=l, A: = |k| = ^, and 

2^2 

< Eq >= iafe^e ^ (0, 0,1) 

= kaijif {0,0,1). (99) 

±af(±^ - “2;, 0)}. (100) 

(e = /r = n = l,k = |k| = —). 

c 


spinning. We also see in (11021) the helicity gradient ra¬ 
dial term, as well as an azimuthal part proportional to 
5{(E* • V)E according to Ea. dMl) . 

No z-component due to 0AM transfer is however ex¬ 
hibited by (11021) . This is a consequence of obviating the 
fundamental conditions: V-E = 0, V-B = 0in the 
paraxial representation (1551) of the Gaussian beam. Of 
course, neither all terms of the decompositions (1551) and 
(1551) are valid for this beam. Hence, although some as¬ 
pects of the force exerted by a TE and TM beam, like 
that of (1551) . are easily described with such representa¬ 
tion of its E and B vectors, other effects of the dynamics 
as discussed in Eg. (11021) . are not adecuately described. 

Gaussian beams satisfying the transversality 
condition 

Paraxial expressions of vortex Gaussian beams that, 
being dz — ikz, satisfy the transversality condition, 
should have = {i/kz)V_L ■ e|^^ [^, e|(^ = 

{Ex \ Ey'^). A first example, that we denote as PGl, has 
the annular spatial field distribution in the AE-plane: 

E«(r) = efc'A'(l,±*, f e±'^[^(l T 1) 

OR d2 

(103) 

Alternatively, a second instance that we denote as PG2, 
is similar to an hypergeometric Gaussian beam in the 
pupil plane: 


As expected, Eos. (1551) - (llOOl) show features similar to 
those seen above for the Bessel beam concerning the in¬ 
fluence of the structure of the beam, in both the longitu¬ 
dinal and tranversal components of these torque parts. 

On the other hand the recoil torques are 

2 J?2 

= = ( 0 , 0 , 1 ); ( 101 ) 


EC)(r) = e|l,±i,(±^ 



(104) 

B(*Vr) = -i-V X eW. 
k 

(105) 

{x,y,z), R=(a;,y,0). 



which are contributed by the electric helicity flow term: 
of Eqs.(15ni) and ([55]). 

When the resultant torque, Eq. dM)) . < r(r) >= C + 
= C + + $( 5 ^ -I- exerted by this beam on the 

particle, is evaluated, one obtains: 

<r(r) >=e2e-^{±^a(“)(0,0,l) 

±ae(±^ - 0)]}- (102) 

Once again, we obtain in the first term of ()102l) the 
spin torque, with its ± sign, along OZ, proportional to 
the absorption cross section, accounting for the particle 


We choose I = 1 and the upper signs in (I103p for the 
PGl beam. According to the first term of (1721) . the spin 
torque < To- > that it induces on the above dielectric 
dipolar particle is: 


< r >,^= e^k^R^e 


1 


1 


27rfc ^ /CzCr^ 
1 


{-y,x,o) 


-(0,0,1)}. (106) 


Whereas by the same token, the PG2 beam gives rise 
to: 


„ _ 2^2 (j(“) 

<r><.= ±ee ^^{(± 



+ (0,0,1)}. (107) 
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R (nm) 


FIG. 3: (Color online). Doped Si sphere of radius 230 nm illu¬ 
minated at A = 1350 nm by either of the Gaussian beams PGl 
[Ea. H103ll ] and PG2 [Ea. dlOlIl ] with LCP transversal com¬ 
ponents. 1 = 1 for PGl, and I = 2 for PG2. cr = 2fim. 

= 1.095 X lO^nm®, al = 6.083 x 10®nm^. Full red: Longi¬ 
tudinal component of < F >2 due to PGl. Full blue: Longitu¬ 
dinal spin torque < Fo- > from PGl, (its magnitude appears 
multiplied by 10®). Broken brown: Longitudinal component 
of < F >2 due to PG2, (its magnitude appears multiplied 
by 14). Broken pink: Longitudinal spin torque < Fo- > from 
PGl, (its magnitude appears multiplied by 10®). The four 
quantities have circular symmetry with an annular spatial 
distribution in the transversal XT-plane, and thus in this 
representation they are even functions if the axis of the radial 
coordinate R is extended to the left of 0 . They are plotted in 
nm® since they are normalized to half the incident intensity 
density: e^. 


For brevity we concentrate on these longitudinal z- 
torques, as this component was the source of the dis¬ 
crepancy encountered with (I88p . They are: 


< r(r) > 2 = e 


Att 



1 )]}( 0 , 0 , 1 ), 


for the PGl beam. And 


9 r 0 "^°^ 

<r(r) >,= ee-T-{± — 
^(;ti±^K}(0,0,i), 


(108) 


(109) 


for the PG2 beam. 

Eos. (11081) and (|109l) contain the topological charge, or 
orbital number, I, coming from the < S > component 
of the orbital momentum Pq, Eq. (jM]) . which describes 
the trajectory of the particle around the vortex following 
the sign of 1. There is no contribution to this I term of 


< r(r) >2 from 3[(E* • V)E] since, for example for PG2, 
2 

the latter yields (±4/cr^)(j/, —x, 0), whose mo¬ 

mentum R X 3[(E* • V)E] contributes to the 
factor of al in (11091) . Of course, again the sign of the inci¬ 
dent spin angular momentum is followed by the spinning 
particle as it absorbes incident energy through . 

Fig. 3 shows the annular spatial distribution of the 
spin and longitudinal total torques on the Si sphere ad¬ 
dressed above, again illuminated at A = 1350 nm by ei¬ 
ther the beam PGl (with Z = 1) or PG2 (with I = 2), 
choosing the upper sign in (Iin3|) and (I104|) . The lon¬ 
gitudinal OX-torque is dominated by the term of 
Eos. (11081) and (11091) . which is several orders of magni¬ 
tude larger than the spin torque < >. The latter 

clearly displays the beam transversal Gaussian shape at 
whose peak the particle is stably placed by the gradient 
force, while it orbits around the vortex and spins on its 
axis. Increasing |Z| enlarges the strength of < P > 2 . 

It is remarkable that while the PG2 beam induces a 
resultant z torque that follows its helicity like the spin 
torque, the interaction of the transversal and longitudinal 
components of the PGl wavefield induces a negative elec¬ 
tromagnetic torque on the particle, which consequently 
orbits around the vortex axis in opposite sense to its spin¬ 
ning. This example, as well as the one above concerning 
a TE Bessel beam, illustrate that in contrast with neg¬ 
ative radiation pressure and tractor beams 3§-^|, the 
production of negative torques is quite common and ex¬ 
ists beyond illumination with special oscillating spatial 
profiles. 

We should remark that in the examples here shown, the 
torque is enhanced by illumination in the electric dipole 
resonance region of this kind of high np particles. By 
contrast, we obtained torques one order of magnitude 
smaller on latex spheres (up = 1.5) of the same size in 
the equivalent region (A = 700 nm). 


CONCLUSIONS 


We have studied the significance and consequences of 
the conservation laws of the spin and orbital angular mo¬ 
menta for the scattering of arbitrary wavefields, satisfy¬ 
ing the transversality condition, by an object, which we 
have illustrated by a magnetodielectric bi-isotropic dipo¬ 
lar particle in the wide sense. We have shown that these 
laws describe the respective contributions of these angu¬ 
lar momenta to the different parts of the torque exerted 
by the field on the particle, which is seen to consist of an 
extinction and a scattering, or recoil, component. The 
latter is shown to cancel the often called intrinsic torque, 
contained in the extinction component and modelled in 
some previous works, without compatibility with energy 
conservation, to account for the optical torque felt by the 
body, and that as shown here plays no role in the resul- 
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tant optical torque, which is actually felt by the particle 
through its absorption cross section. This latter quantity 
is what remains from the recoil torque component after 
the above mentioned cancellation. 

Each of these two laws describe a half of the recoil 
torque acting on the object, The transfer of spin and or¬ 
bital angular momenta, on the other hand, arises while 
they are extinguished from the incident field. This is de¬ 
scribed by the torque extinction parts due to the interfer¬ 
ence of the incident and scattered fields. If the incident 
wave is plane, and thus it has no orbital angular momen¬ 
tum, its transfer characterized by the recoil torque is a 
manifestation of the spin-orbit interaction. 

The role of the spatial structure of the incident field in 
the torque has been established, showing that the con¬ 
tribution from the spin angular momentum conservation 
is twice that from the orbital angular momentum con¬ 
tinuity law. While the spinning of the particle through 
absorption of the incident energy always follows the in¬ 
cident helicity, its orbiting may result opposite to the 
incident spin, this giving rise to a resultant negative op¬ 
tical torque on a single particle. I have illustrated this 
with Bessel and Gaussian beams. 

The electromagnetic torque admits a decomposition 
into conservative and non-conservative components in 
which the helicity and its flow play a role analogous to 
that of the energy and the Poynting vector in the optical 
force. In particular, the gradient torque has the potential 
to add a new variable to optical tweezer set-ups. 

In addition to contributing to an understanding of the 
physical mechanisms involved in the transfer of angular 
momentum from the field to a body, ruled by these con¬ 
servation laws, and their several consequences, this study 
opens new avenues of research both from the experimen¬ 
tal and theoretical points of view that, with the improve¬ 
ment of particle manipulation techniques, adds new de¬ 
grees of freedom in the analysisis and control of possible 
fundamental phhenomena and their applications. 
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